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Abstract. - We study the dynamics of a single polymer subject to thermal fluctuations in 
a linear shear flow. The polymer is modeled as a finitely extendable nonlinear elastic (FENE) 
dumbbell. Both orientation and elongation dynamics are investigated numerically as a function 
of the shear strength, by means of a new efficient integration algorithm. The results are in 
agreement with recent experiments. 



Introduction. Nowadays, thanks to the development of effective experimental tech- 
niques, it is possible to follow the motion of a single macromolecule in a flow, either laminar 
or turbulent [1-8,11-18]. This is of crucial importance for applications in polymer process- 
ing [19] and biophysics [1]. Dynamical properties of biomolecules have been explored in detail 
(see e.g. [1-5,11-14] for DNA and [15] for chromatin) and protein-macromolecule interactions 
have been studied [16-18]. 

The formulation of theoretical models (see e.g. [19]) able to reproduce qualitatively and quan- 
titatively these measurements represents an important step towards the understanding of 
single-molecule biophysics. An extensive analysis of single polymer dynamics in simple flows 
has been conducted in a series of papers by Chu and coworkers, Larson and coworkers and 
Shaqfch and coworkers (see [1-8] and references therein). Here we mention in particular two 
recent papers where the statistics of orientation and conformation of long-chain molecules in 
linear shear flows has been studied in great detail, with a direct comparison with numerical 
models [9, 10]. An intrinsic difficulty is represented by the large number of degrees of freedom 
required to describe the polymer conformation, and thus its dynamics. Nonetheless, non- 
trivial aspects of polymer-fluid interactions may be accounted for and even explained at a 
semi-quantitative level by means of simple, few degrees of freedom models. One of the sim- 
plest, yet reliable, model is the finitely extendable nonlinear elastic dumbbell (FENE) [19]. 
The polymer is described by its end-to-end distance vector R and the microphysical properties 
are essentially dumped into two parameters: I/7, the longest elastic relaxation time of the 
macromolecule, and £, its friction coefficient with the surrounding solvent. Nonlinear clastic 
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Fig. 1 - A sketch of the geometry of polymer motion in a linear shear flow. 

effects must be accounted for whenever the polymer is considerably stretched, as is the case of 
shear flows [4]. The geometry of this problem is depicted in Fig.QJ The polymer spends a large 
fraction of time in elongated configurations along the shear direction. In the following we will 
present numerical results about end-to-end orientation, elongation and about the statistics 
of tumbling times. The latter is defined as the time spent between two successive "flips" of 
the polymer ends (see Fig. Tumbling can occur via different pathways, e.g. passing by 
a coiled state or through folded configurations: those details cannot be addressed within the 
single-dumbbell model and will be the subject of future study. 

From experiments to numerics. - Single polymer orientation and tumbling dynamics was 
recently studied experimentally by Steinberg et al. [12] with a 10 -3 ppm solution of A— DNA 
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Fig. 2 - Four possible stages of a tumbling event. 
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molecules labeled with fluorescence methods (but see also the experiments presented in [9,10]). 
To resolve the angular dynamics two different flow configurations were used: one was generated 
by two discs one of which was rotating with uniform angular velocity Q and the other flow was 
obtained by a boundary layer in a micro-channel produced with the soft lithography method 
(see [12] and references therein for details). 

To reproduce the physical situation of [12] we studied a FENE dumbbell in a simple shear flow 
v = (sy, 0, 0), sensing thermal fluctuations (see Fig.^l. The equation describing the evolution 
of the end-to-end vector of the polymer is: 



where R = KgT/H , 7 = 4iJ/£, rj is a three-dimensional white noise with zero mean and 
correlation (ViityVj (t')) — ~ t'), R m is the maximum length of the polymer, Kg is the 

Boltzmann constant, T is the temperature, £ is the isotropic drag coefficient, and H the spring 
constant. 

Even if the single-FENE-dumbbell model does not reproduce precisely the behavior of real 
molecules [8], we can set the parameters of our model, Rq, R m , 7, as close as possible to 
their corresponding experimental values [4,12]: we choose Rq ~ 1 fim,, R m /Ro ~ 21, 7 ~ 
0.01s- 1 ~ Is" 1 . 

The orientation dynamics has been investigated for rigid spheroid by Hinch and Leal [20] . As 
for polymers, at large Weissenberg numbers Wi = s/7, where s is the shear rate, the basic 
ingredients of the polymer dynamics can be summarized as follows [7,20-25]: due to the 
shear flow the polymer tends to reach the unstable equilibrium configuration where it is fully 
extended along the shear direction. In polar coordinates (R, 9, 0) = (Rm, 0, or %). The effect 
of thermal noise is to drive the polymer away from this configuration. The most probable value 
of 9 is zero, due to the symmetry of the dynamics along the z axis. However, large fluctuations 
in the off-shear-plane angle can occur. The most probable value for <f> will be slightly larger 
than 0, or 7r: the symmetry-breaking effect of shear causes the polymers to "hesitate" for 
some time before crossing the x axis and then give rise to a tumbling event. Few results can 
be obtained analytically for this model, except for the linear case where R m /Ro — > 00. 



Numerical algorithm. - Several numerical methods have been proposed to simulate poly- 
mer dynamics (see for example [26]). A commonly encountered problem with nonlinear elastic 
models is the loss of accuracy close to the singularity R — » R m . In order to overcome this 
problem it is possible to perform a change of variables in the vicinity of R m that removes the 
singularity and allows to use a straightforward time-marching scheme. This method can be 
easily extended to other nonlinear models [27] as well as to other flows. 

Eq. can be solved by any stochastic discretization scheme (Euler-Ito in our case) in the 
region R < Rthr, where Rthr is a fraction of R m , say 0.5i? m . Whenever R exceeds the thresh- 
old we switch to polar variables (R, n) , where h is the unity vector describing the orientation 
of the polymer hi = Ri/R, and then to the new variables (z, h), where 
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Fig. 3 - The PDF of elongation plotted for Wi = 0, 1, 5, 10, 40, 200 (from A to F). 



This relation can be easily inverted to give R as a function of z. After computing all the 
contact terms in the Ito convention we have the following equations for (z,n): 
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which is regular in the neighborhood of R = R m , i.e. z = 0. 



Results. - The Probability Density Function (PDF) of the modulus of the conformation 
vector depends strongly on Wi. At sufficiently small Wi sa 1 the statistics does not differ much 
from the linear elastic case, since R <C R m . The PDF in the FENE case can be computed 
analytically only in asymptotic regimes [25,28]. The numerical result is shown in Fig. [3] for 
several Weissenberg numbers. At very large Wi the elongation PDF presents a peak with 
height scaling as Wi 2 ^ 3 and width as Wi~ 2//3 (not shown). This result is in agreement with 
the predictions of Ref. [23] . As a side remark we notice that experimental measurements of 
the elongation PDF at Wi as large as Wi = 76 do not display a peak near R m (see [8] and 
fig. 5 of [4]) as well as numerical results of multi-beads models (see figs. 4, 5 in [8] and the 
discussion therein). 

The orientation of polymers follows the qualitative picture drawn in the linear elastic 
case [24,25], even though there appear important quantitative differences. The PDF of the 
in-shear-plane angle (fi is shown in Fig. 0] The probability is concentrated in the vicinity of 
(j> = 0, 7r with a peak width at half height (fit, whose dependency on Wi is shown in Fig. EI The 
case of a linear elastic dumbbell R m = oo is shown for comparison. The angle <fi t decreases 
with Wi in both cases, i.e. the larger is Wi the narrower is the region around the x axis where 
the polymer spends most of its time. The scaling can be derived by simple physical arguments: 
following Chertkov et al. [22] the evolution equation for <fi in the region (fi <C 1 is approxi- 
mated by by d t (fi = —s(fi 2 + y/jR^/R 2 r/^, where rj^ is a white noise. Thus (fit can be estimated 
balancing shear and noise terms in the right hand side terms, i.e. (fit ~ Wi~ 1 ^ 3 (i?o/^) 2 ^ 3 - 
At large Wi, for a linear elastic dumbbell one has R oc Wi, yielding (fit ~ Wi -1 , whereas for 
a nonlinear elastic force one estimates R ~ R m to find <fi t ~ Wi -1 / 3 . The tails of the PDF 
follow closely the distribution sin -2 <fi dictated by the shear (see Fig. 0J. 
The agreement with the experiments is very good [12]: the scaling in the tail follows sin~ 2 (fi, 



A. Celani et al: Polymers in linear shear flow 



5 




Fig. 4 Fig. 5 

Fig. 4 - The PDF of the angle cf> plotted with sin -2 0. In the inset the PDF in linear scale plotted 
for Wi = 1,5,40. 

Fig. 5 - The behavior of <f>t as a function of Wi. 

and the dependence of </>t on Wi is close to the theoretical prediction already for Wi = 25. 

The marginal PDF of the angle 9 is presented in Fig. The tails decay as 9~ 2 , with 
a scaling range increasing with Wi. The algebraic behavior has been observed in [12] for 
Wi = 17.6, and even if Wi is not very high the agreement is remarkable. The probability 
density of 9 for small angles cj) ~ 4> t , or equivalently the joint PDF P(9,cf> = 0), shows a 
neat power law close to 9~ 3 for 9 3> 9 t . This nontrivial scaling behaviour has been predicted 
theoretically and observed numerically for the linear elastic case in Refs. [22,24,25]. The 
width of the peak of the P{9) at half height, 9t, decays as Wi -1 ^ 3 for the nonlinear elastic 
case. The agreement with experimental data is perfect [12]. 

Note that the cross-over between the linear elastic case 6 t ~ t ~ Wr 1 and the FENE case 

1/3 

<t>t ~ 9 t r~j Wi ' occur at Wi ~ 1, as measured experimentally. 

For what concerns the tumbling times statistics there are two possible definitions [25]: (i) 
given an appropriate threshold value in R that defines the coiled state for the polymer, one 
can compute the time spent during two successive coiled states [12]. (ii) one can consider the 
time between two subsequent crossings of the plane 4> = n/2. 

Both definitions are ambiguous for small values of Wi, i.e. when the polymer spends most 
of its time in a coiled state. For large tumbling times the PDF is exponential for both 
definitions of r, P(r) ~ exp [—-Br] (see Fig.|HJ). This is a robust feature of this phenomenon [12, 
22]. Experimental measurements of the tumbling time are possible only following the first 
definition, due to lack of angular resolution [12]. 

The exponent E of the tail in the linear elastic case is inversely proportional to the relaxation 
time of the polymer 7 -1 . In the FENE case there is a non trivial dependence on Wi, as shown 
in Fig. El The scaling of the typical tumbling time r t ~ E^ 1 can be estimated at large Wi as 
follows: the angular motion in the region <p ~ <p t is driven by the thermal noise and is therefore 
diffusive. The diffusion coefficient is D = j(Ro/R) 2 and therefore r t ~ 4> 2 /D. Substituting 
R oc i?oWi for the linear spring model and R ~ R m for the FENE model one obtains for E/j 
the scalings Wi° and Wi 2 / 3 , respectively. 

The behavior of the PDF at r < r f is model dependent and should not be considered as 
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Fig. 6 Fig. 7 

Fig. 6 - The PDF of the angle 9 plotted against 9~ 2 . In the inset the PDF in linear scale for 
Wi = 20, 40, 100. 

Fig. 7 - The behavior of 9 t as a function of Wi. 



relevant (see e.g. [25]). In experiments [12] the exponential tail of the PDF is observed and 
the dependence of r t on Wi is in accordance with theoretical arguments and numerical results. 

Conclusions. - We studied the dynamics of a single FENE polymer immersed in a simple 
shear flow with thermal noise. The statistics of orientation, elongation and tumbling of the 
polymer have been analyzed in comparison with experimental measurements [12], previous 
numerical simulations [8,25], and theoretical expectations [22-24]. Even if the large variety of 
conformations of real polymer can not be explored within such a simple model, single-FENE- 
dumbbcll can reproduce semiquantitatively several aspects of the behavior of real polymers. 




Fig. 8 - The PDF of tumbling times: tr is the time elapsed between two nonadjacent coiled states, 
defined as the states where R is smaller than XRo (A = 1.5 in this figure); t$ is defined as the time 
between two rotations of 7r in the angle 4>. Here, Wi = 76. 

Fig. 9 - The exponent E rescaled with the relaxation time as a function of Wi in the linear case and 
in the FENE case. 
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